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We evaluate the charge and longitudinal spin response functions of a two-dimensional electron gas 
("^ ' with e^/r interactions in an arbitrary state of spin polarization, using a structurally self-consistent 

^3 ' approach to treat exchange and correlations. From the results we assess the nature of the magnetic 

Cn , order in the electronic ground state in zero magnetic field as a function of electron density. We find 

that states of partial spin polarization are thermodynamically unstable at all values of the coupling 
strength and that a first-order phase transition occurs with increasing coupling strength from the 
magnetically disorderd (paramagnetic) phase to the fully spin-polarized (ferromagnetic) phase. This 
behavior is in qualitative agreement with diffusion Monte Carlo data, although the location of the 
phase transition is underestimated in our calculations. 
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I. INTRODUCTION 



''O ■ Electronic systems that exhibit quasi-two-dimensional (2D) behavior have attracted great interest for a number of 

b^ \ years. They include electrons in semiconductor heterojunctions and quantum wells and in inversion and accumulation 

Q ■ layers at semiconductor-insulator interfaces^. The basic theoretical model for such systems is the gas of electrons 

(EG) interacting with the e^/r law and moving in a plane over a uniform neutralizing background of positive charge. 

'""' i It is a general feature of interacting electron fluids in both 2D and 3D that spin fluctuations and their correlations 

JC , become more important as the coupling strength increases with decreasing electron density. This fact has been most 

T-H strikingly emphasized by quantum simulation work. In the 3D fluid these studiesoO have revealed a continuous 

Jl'J transition from the paramagnetic to the ferromagnetic state, before a first-order transition to a ferromagnetic Wigner 

f~^ ■ crystal occurs. In the 2D EG, instead, similar studiesErEl have found no evidence for partially spin-polarized ground 

states, but have indicated that a first-order transition occurs from the paramagnetic to the ferromagnetic fluid before 

crystallization. 

C^ ' A theoretical study of the magnetic phase behavior of the 2D electron fluid is the purpose of the present work. 

^ ■ Previous calculations of charge and spin correlations in the 2D EGQ have referred to the paramagnetic state only. We 
i-rt consider instead, at a chosen set of values for the Coulomb coupling strength, a 2D EG in an arbitrarily chosen state 
^ , of spin polarization and evaluate the whole matrix of dynamic susceptibilities for its charge and longitudinal spin 
r^ I response. From this viewpoint we can assess both the absolute dynamic stability and the relative thermodynamic 
stability of the various states of spin polarization. The basic approximation underlying our calculations is an immediate 
extension of the so-called STLS approachQ to treat exchange and correlations through self-consistency between linear 
susceptibilities and pair distribution functions. This is a necessary starting point to the problem at hand and allows, 
Cd [ as we shall see, qualitative understanding to be gained on the magnetic phase behavior of the 2D EG. 

The contents of the paper are briefly as follows. Section O presents the essential definitions of linear response 
functions and their expressions through local field factors allowing in an unrestricted way for exchange and short- 
range correlations between the three types of electron pairs in a spin-polarized stateH. The approximate expressions 
of the local field factors embodying self-consistency with the liquid-structure pair functions are also given there. 



Section |II| reports the main results of the self-consistent calculation of structure and response and then focuses 
on the numerical results for the ground state energy as a function of electron density and spin polarization. We 
show that, starting from the usual paramagnetic ground state in the weak-coupling regime where we find very good 



agreement with the accurate results reported by Yarlagadda and Giuhanfl for the spin susceptibihty, the states of 
partial spin polarization are squeezed out before a crossing occurs between the energy curves of the paramagnetic 
and ferromagnetic states. The transition between these two states is found to occur, however, at a lower value of the 



coupling strength than indicated by the quantum simulations. Finally Section IV summarizes our main conclusions 
and gives suggestions for further work. 

II. THEORETICAL APPROACH 

We consider a 2D EG modelling a system of electronic carriers with effective mass m* in a semiconductor with 
background dielectric constant eo and interacting via the Coulomb potential Vq — 2T:e'^ /e^q. The thermodynamic 
state parameters of the EG at zero temperature are taken as the areal densities n^ and n^ of spin-up and spin-down 
electrons or equivalently the total density n = n| + nj^ and the spin polarization C, = (n| — ni)/n. The coupling 
strength parameter is r^ = (7rna|j^)~^/^, where a*^ = h^eo/m*e'^ is the effective Bohr radius. 

The charge and longitudinal spin response of the EG to an external electromagnetic field is described in the linear 
regime by a set of susceptibilities Xctct'(9, w), the subscripts a and a' being spin indices (t or ['p. These arc defined 
through the double Fourier transform of the linear response functions of the EG to a position and time dependent 
electromagnetic field, which are given by 

K,„,{\v - r'l, t - t') - -ine{t - t') ([p,(r, t), p,,(r', t')]) . (1) 

Here, 6{t) is the Heaviside step function and p„{r,t) is the spin density operator in the Hcisenberg representation. 
Linear combination of the Xcja' susceptibilities give the charge susceptibility Xcc and the spin susceptibility Xss as 
well as the off-diagonal susceptibilities Xcs and Xsc- Specifically, leaving implicit the dependence on q and w, we have 
Xcc = e^Y^a^'Xaa'-, Xss = "7^ Eaa' sgn(cro-')Xo-o-' , Xcs = 67 E^^' sgn(cr')x<T<T' and Xsc = e7Eaa'Sgn(cr)x<T<T', with 
7 — 9[J^bI'2 where g is the Lande factor and hb the Bohr magneton. 

The partial structure factors Saa' (?) are related to the linear susceptibilities by the fluctuation-dissipation theorem, 
which can be written in the form 

Saa'{q)^ , / Xaa'iq^iu). (2) 

The integral in Eq. (y) is taken along the imaginary frequency axis, the contribution from the plasmon pole being 
easily handled numerically by this method (see de Freitas et aln). The pair distribution functions ga-a-'{r) are obtained 
from Eq. (0) by Fourier transform. 
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gaa' (r) = 1 + ^ / 7^—2 [Saa' {q) - Saa'] exp(iq • r) (3) 



and the ground state energy E(rs, C) as a function of electron density and spin polarization can in turn be calculated 
from these structural functions. The results is 



E{r,,0=^-^ + ^^^^^ I dr'Mr'sX), (4) 



in units of the effective Rydberg Ryd* — e^/2eoa*B- The first term in Eq. (Bl) is the free-electron kinetic energy and 
the second term is the exchange-correlation energy E^c, with 7 defined by 

/•OO 

7(r„C)= / dq[S,,{q)~l] (5) 



and Scc{q) = [^T'^'ttI?) + ^iSuil) + 2(n|n|)^'^S'||(g)]/n being the charge-charge partial structure factor. One can 
extract the correlation energy E^ from E^c by subtracting the Hartree-Fock exchange part. We have 
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(6) 



The local field factors Ga-a' [l, w) are introduced by writing an effective one-electron Hamiltonian in which an 
electron with spin a experiences effective fields embodying exchange and correlation with the surrounding EGS. In 
the following we adopt an STLS-like approximation, in which the bare electron-electron interaction Vq is replaced by 
an effective spin-dependent static interaction ■^^.^'(g) = Vq[l — Go- o-' ((?)]. The corresponding expression for the linear 
susceptibilities is 



Xaa'{q,Uj) 



XO<T(g,w) [S^^' + V<7a''ipaS'{q)X0S'iQ,>^)] 



A(g,w) 



(7) 



where a denotes the spin orientation opposite to a, 77^0-' 
the ideal 2D Fermi gas and A{q,uj) is given by 



(_1)'5<T^'^ Xoct('Z, ijj) are the (diagonal) susceptibilities of 



A{q,uj) = [i-ipf^{q)xoiiq,uj)] [1 - V'ii('7)xoi('?, t^)] ^ ■iPtii<i)'^ui9)xm{q,^)xoi{q,^) ■ 



(8) 



The expressions for the real and imaginary part of the ideal-gas susceptibility arc well known (see e.g. Ref. [ H). 
From these the expression for xoaiqiiu) on the imaginary frequency axis is obtained for use in Eq. (m) as 
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where we have defined 
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with kpa = (47rncr)^/^. Finally, the local field factors are determined from the partial structure factors through the 
closure relationfl 
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S^„,{\k~q\)-5„„,]. (11) 

I, leading at convergence to numerical results 



^riaUa' J (27r)2 k 

A self-consistency cycle can thus be set up by using Eq. (|^) in Eq, 
for the partial structure factors and for the local field factors. 

With the ground state energy as a function of r^ and C at our disposal we can calculate the compressibility and 
the spin susceptibility of the EG. The compressibility k is given by the partial derivatives of the ground state energy 
with respect to the density: 
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where we have used the definition of the pressure P = ~{nra/2)dE /drs\Q. The spin susceptibility Xs = liniq^o Xss{(ii 0) 
can similary be obtained from the second derivative of the ground state energy with respect to magnetization: 

XP m* d'^E 

Xs 



in' oe 



(13) 



The quantity xp ~ rn* g'^ /j^^ / Anh is the Pauli spin susceptibility of the ideal 2D Fermi gas. An alternative route to 
the compressibility and the spin susceptibility is through the long-wavelenght limit of the local field factorsQ, as will 
be discussed for the paramagnetic state in Section III. B below. 

The function A (9,0;) also carries valuable information about the collective excitations and the dynamical stability 
of the system. The plasmon dispersion relation LOp{q) is found by solving the equation A(g,tj) = 0. The existence of 



a solution for the equation A{q, w) = when w is purely imaginary is a signature of a dynamical instability in the 
system. We shall meet this instability in Section III. C. 

As a final remark, the matrix of charge-spin susceptibilities becomes diagonal in the paramagnetic state, yielding 



where xo(9i<^) = 2xot(9,^) = ^^Xoiil^^)- Equations (|l3) and (0) can be rewritten as 
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and 
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where kq = ^rrf /2 is the ideal-gas compressibility in units of a*g / Ryd* . 

III. NUMERICAL RESULTS 

In this section we report our numerical results for the structure and local field functions and for the ground-state 
energy over relevant ranges of the coupling-strength parameter r^ and of the magnetization fraction (. We also report 
the plasmon dispersion relation and discuss the stability of the 2D EG as a function of (. 

A. CORRELATION FUNCTIONS AND LOCAL FIELDS 



Figures and show the behavior of the partial structure factors as functions of q/kp] , first with varying C, at 
rs = 2 and then with varying r^ at C = 0.5. Starting from the left in each figure, the three sets of curves give Sn^q), 
Siiiq) and S^iiq) = Si^{q). 

It is seen from Figure |l| that with increasing ^ away from the paramagnetic state {i.e. with the induction of 
spin polarization in the up direction) the function S-\-\{q/kp-^) remains essentially the same while Sii{q/kF-\) shrinks 
towards the vertical axis. These behaviors descend from the fact that these structural functions mainly reflect the 
Pauli exchange hole between parallel-spin electrons. Accordingly , these structure factors are almost independent of 
the Coulomb coupling strength, as is seen from Figure 0. The cross structure factor S'ji{q/kp^), on the other hand, 
shows a trough due to correlations from Coulomb repulsions, which deepens with increasing coupling strength in 
Figure g. The shift of the trough of S'ii{q/kF-\) in Figure |l|, as well as the shrinking of Sii{q/kp'i), are mainly due to 
the increase of the scaling wave numbers kp^ with increasing ^ at constant total electron density. 

The behavior of the Pauli exchange and Coulomb holes is directly seen in the pair distribution functions ^ctct' 
and is illustrated with varying C at I's = 2 in Figure and with varying r^ at C — 0.5 in Figure 0. The holes in 
5ll(r) and gii{r) show little dependence on both coupling strength and spin polarization (except for scaling of the 
unit of length). They do not exactly vanish as rkp^ — > 0, which evidently is a quantitative defect of the present 
approximate approach. The Coulomb hole in g]i{r) — .9j,t(?'), on the other hand, shows appreciable dependence on 
coupling strength (as expected) and on spin polarization. Except for the region of small rkp^, an increase in spin 
polarization broadens the Coulomb hole and thus brings about an effective decrease in the Coulomb repulsive energy. 
This ultimately leads to stabilization of the ferromagnetic state. The very weak dependence of .9tx(0) on the state of 
spin polarization has already been discussed elsewherdla. 

Finally, Figure |^ shows that the local field factors Gaa' (q) have a rather weak dependence on spin polarization 
at fixed r^ = 2. The asymptotic values of these functions for q — > oo are determined by the relations 6*11(00) — 



G'xx(oo) = 1 and G||(oo) = 1 — (7|j(0)eI: these relations are rather well satisfied in our approximate approach. The 
slopes of Gctct' (q) in the long wavelength limit are instead related to the compressibility and spin susceptibility and 
will be discussed below. 

B. GROUND STATE ENERGY 

We evaluate the ground state energy E as a function of rg and C from Eq. (Q) . Figure y reports £■ as a function of r^ 
for various values of C in the range to 1. The paramagnetic state has the lowest ideal-gas kinetic energy and remains 
thermodynamically stable over the whole range of coupling strength illustrated in Figure o, up to r^ = 4. The energy 
difference from the ferromagnetic state is steadily decreasing and the energy curves corresponding to intermediate 
values of C are being squeezed together above their various minima. In fact, as we shall show in full detail in the 
next section, the states with partial spin polarization become dynamically unstable in the upper part of this range of 
Tg. Only the paramagnetic and fully spin polarized states remain dynamically stable above Ts = 4 and their energy 
curves cross at Tg = 5.5, as is seen from Figure R. 

The present approximate approach thus predicts a first-order quantum phase transition occurring from the para- 
magnetic to the ferromagnetic state at r^ = 5.5 in the 2D EG at zero temperature. Above their crossing, the energy 
curves of the two states remain very close to each other over the whole range of r^ illustrated in Figure 0, up to 
Ts = 10. There is, therefore, a large uncertainty attached to our prediction of the location of the phase transition. 

Figure g reports a number of different results for the spin susceptibility Xs of the 2D EG in the paramagnetic state, 
in units of the ideal Pauli susceptibility xp, as a function of r^ up to r^ = 4. The values of Xs/xp that we obtain by 
differentiation of the ground state energy from Eq. (Oh (denoted as STLS2 in Figure ph are in reasonable agreement 
with the values extracted from the available quantum Monte Carlo datati, which are shown as crosses in Figure & 
We also obtain close agreement with the theoretical results of Yarlagadda and GiulianiB at low values of rs , although 
these values of Xs/xp increase much too rapidly with coupling strength beyond the range of r^ studied in the original 
work of these authors. A general comment applying to both their results and ours is that the role of exchange is 
somewhat overemphasized relative to the Monte Carlo data. The location of the phase transition is accordingly placed 
at a lower Vs , as we have already seen in Figure M. 

There is, however, a large violation of the susceptibility sum rule in our approach: the value of Xs/xp that we 
obtain from the slopes of the local field factors at long wavelengths according to 

^ - 1 - lim[GTT(g) +Gu(g) - Guiq) - G^^iq)]/qa*s (17) 

Xs 1^0 

(denoted as STLSl in Figure 0) are vastly different from those obtained from the ground state energy according 
to Eq. (O) and indeed their dependence on Ts is much too close to that following from the simple Hartree-Fock 
approximation. 

Finally, the results that we obtain for the compressibility of the paramagnetic state by density differentiation of 
the ground state energy according to Eq. (|lj) are in excellent quantitative agreement with the Monte Carlo dataD'El. 
There is again a large violation of the compressibility sum rule, this being a well known defect of the STLS approachQ. 

C. PLASMON DISPERSION AND DYNAMIC INSTABILITY OF PARTIALLY POLARIZED STATES 

The plasmon excitation energy ujp{q) is evaluated from the root of the equation l^{q^Ljp{q)) = 0, A((7,a;) being as 
given in Eq. (@) . The results are shown in Figure at Ts = 2 for various values of the spin polarization, superposed 
on the particle-hole continuum for the majority spin population. It can be shown analytically from Eq. (S) that the 
plasmon dispersion relation at long wavelength is given by u}p{q) —f ^j2Ti{n^ + nije^q/ t^m* . 



We turn next to examine the behavior of the function A{q,u!) on the imaginary axis w — iu. We notice first that 
the energy change AE associated with a small modulation Sua- (q) of the spin populations away from a homogeneous 
state with densities n| and n^ is given by 

AE=-Y,J^X-^,{q)SnAcDSn,'ici). (18) 

(T.(j' 

This quantity must be positive if the undeformed spin state is stable and, using Eq. (^ and the fact that Xot{QjO) is 
negative for all values of q, we find that the stability condition can be written as 

A(g,0)>0. (19) 

In fact, it can be shown analytically from Eq. (||) that in the paramagnetic state the stability condition in Eq. ( [l9| ) 
reduces at long wavelengths to the condition that the spin susceptibility defined in Eq. (|l7| ) be positive. 

Figure nO reports the function A{q, iu) plotted against u a.t C, = 0.5 and various values of Vg for a small value of q 
{q — O.lkpi). It is evident that this state of partial spin polarization is becoming dynamically unstable at r^ « 3, as 
signalled by a change in sign of the function A(q, iu) at low wave number and frequency. 

In view of the violation of the susceptibility sum rule in the present theory as discussed in Section III. B wc should, 
however, return to Figures o and ^ and assess the thermodynamic stability of the spin-polarization states for the 2D 
EG by the well tested and reliable results reported these. That is, states of partial spin polarization are never stable 
and the system goes from paramagnetic to ferromagnetic by a first-order phase transition. 

As a final remark, we find no dynamic instability in the charge-charge response function of the 2D EG in cither the 
paramagnetic or the ferromagnetic state over a wide range of Vs (up to r^. = 30). 

IV. SUMMARY AND FUTURE DIRECTIONS 

In summary, we have evaluated the linear charge and longitudinal spin susceptibilities of a 2D electron gas in 
various states of spin polarization over a range of electron densities, using a structurally selfconsistent approach to 
treat exchange and correlations. The main results of the calculations are the structural pair functions and hence the 
energy of the system as a function of electron density and spin polarization. These results have in particular been 
tested for the paramagnetic state against previous theoretical results and quantum Monte Carlo simulations of the 
spin susceptibility and of the compressibility. 

Prom the calculated energy function we have found that the 2D electron gas undergoes a first-order phase tran- 
sition from the paramagnetic to the ferromagnetic state with increasing coupling strength, without passing through 
intermediate states of partial magnetization. This result agrees with the scenario coming from quantum Monte Carlo 
simulations, although our predicted location of the phase transition appears to be underestimated. 

We have stressed that the main defect of our approach is the lack of thermodynamic consistency in the compress- 
ibility and spin susceptibility, i.e. the disagreement between the values obtained by differentiation from the energy 
function and those derived from the appropriate limit of the microscopic susceptibilities. The latter show dynamic 
instabilities emerging in the states of partial magnetization as the coupling strength increases. These inconsistencies 
point the way to possible improvements in the theoretical approach through a refinement of the closure relation in 
Eq. (|ll|). As is suggested by early work on the 3D paramagnetic electron gaaZHlil, thermodynamic consistency may 
be imposed by allowing for the density and spin-polarization dependence of the partial structure factors used in the 
evaluation of the local fields in the microscopic susceptibilities. 

ACKNOWLEDGMENTS 

This work was partially supported by MURST through the PRIN1999 program. We thank Dr. M. Polini for useful 
discussions. 



^ T. Ando, A. B. Fowler and F. Stern, Rev. Mod. Phys. 54, 437 (1982). 

2 D. M. Ceperly and B. J. Alder, Phys. Rev. Lett. 45, 566 (1980); B. J. Alder, D. M. Ceperly and E. L. Pollock, Int. J. Quant. 

Chem. 16, 49 (1982). 
^ G. Ortiz, M. Harris and P. Ballone, Phys. Rev. Lett. 82, 5317 (1999). 
■* B.Tanatar and D. M. Ceperly, Phys. Rev. B 39, 5005 (1989); S. Moroni, D. M. Ceperly and G. Senatore, Phys. Rev. Lett. 

69, 1837 (1992). 
^ F. Rapisarda and G. Senatore, Austr. J. Phys. 49, 161 (1996); D. Varsano, S. Moroni and G. Senatore, Europhys. Lett. 53, 

348 (2001). 
''M. Jonson, J. Phys. C 9, 3055 (1976); M. Rajagopal, Phys. Rev. B 15, 4264 (1977); U. de Freitas, L. C. loriatti and N. 

Studart, J. Phys. C 20, 5983 (1987); D. Neilson, L. Swierkowski, A. Sjolander and J. Szymanski, Phys. Rev. B 44, 6291 

(1991); A. Gold and L. Calmels, Phys. Rev. B 48, 11622 (1993); Z. Lenac and M. Sunjic, Phys. Rev. B 50, 10792 (1994); 

R. K. Moudgil, P. K. Ahluwaha and K. N. Pathak, Phys. Rev. B 51, 1575 (1995) and 52, 11945 (1995); B. Tanatar and N. 

Mutluay, Eur. Phys. J. B 1, 409 (1998). 
'' K. S. Singwi, M. P. Tosi, R. H. Land and A. Sjolander, Phys. Rev. 176, 589 (1968). 

® D. C. Marinescu and J. J. Quinn, Phys. Rev. B 56, 1114 (1997); M. Polini and M. P. Tosi, Phys. Rev. B 63, 045118 (2001). 
^ S. Yarlagadda and G. F. Giuliani, Phys. Rev. B 40, 5432 (1988). 

^° M. Polini, G. Sica, B. Davoudi and M. P. Tosi, J. Phys.: Condens. Matter 13, 3591 (2001). 
" P. Vashishta and K. S. Singwi, Phys. Rev. B 6, 875 (1972). 



g 04 




FIG. 1. Static structure factors as functions oiq/kp^ airs — 2 and (" = 0,0.4and0.8. From left to right: 5'|x(g), S^-\(q) and 
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FIG. 2. Static structure factors as functions of g/fc^t at ^ = 0.5 and r^ = land 3. From left to right: S'||(g), S^^{q) and 




FIG. 3. Pair distribution functions versus rkp^ at rs = 2 and i^ = 0, 0.2, andO.4. From left to right: g-jiir), <7tt(»") and 
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FIG. 4. Pair distribution functions versus rkp] &t C, = 0.5 and Vs = 1 and2. From left to right: g-\[{r), g]-\{r) and gu(r) 
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FIG. 5. Local field factors Gaa-'il) as function oi q/kp^ at rs = 2 and (^ — 0, 0.2 and 0.4. The meaning of the curves is shown 
in the inset. 
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FIG. 6. Ground state energy as a function of rs at (^ — 0, 0.3, 0.5, 0.7 and 1. 
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FIG. 7. Ground state energy as a function of Vs for the paramagnetic and ferromagnetic states, ^ — Oand 1 respectively. 
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FIG. 8. Spin susceptibility Xs/xp in the paramagnetic state as a function of r^ from the Hartree-Fock approximation (HF), 
the Yarlagadda-Giuhani calculation (YG), the present calculation according to Eq. (fo) (STLSl) and to Eq.(Fl3) (STLS2), and 
from Monte Carlo data (crosses) 
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FIG. 9. Plasmon excitation energy as a function of q/kp^ at r^ = 2 and C, =Q, 0.2, 0.4, 0.6 and 0.8 (from top to bottom). 
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FIG. 10. A{q,iu) as a function of u &t i^ — 0.5 and rs — 1, 2, 3 and 4. 
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